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Abstract 

Symmetry breaking surface fields give rise to nontrivial and long-ranged 
order parameter profiles for critical systems such as fluids, alloys or magnets 
confined to wedges. We discuss the properties of the corresponding universal 
scaling functions of the order parameter profile and the two-point correlation 
function and determine the critical exponents r}\\ and rj± for the so-called 
normal transition. 
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I. INTRODUCTION 



Advanced experimental techniques have emerged which allow one to endow solid sur- 
faces with stable geometrical structures which display a well-defined design on the scale 
of nanometers or micrometers. Lateral geometric structures can be formed by using var- 
ious lithographic techniques such as, e.g., holographic X-ray (LIGA) 0, soft 0, and 
nanosphere lithography . These microfabrication techniques provide routes to high-quality 
patterns and structures with lateral dimensions down to tens of nm. These structures are 
either periodic in one lateral direction consisting of grooves with various shapes of the cross 
section (e.g., wedgelike) or they display periodicity in both lateral dimensions. 

These man-made surfaces offer a wide range of possible applications if they are exposed 
to fluids. In that case the surfaces act as a template with designed topography which 
imposes specified lateral structures on a fiuid. For example in the context of microfiuidics 
I^J^ these geometrical structures can be used as guiding systems in order to deliver tiny 
amounts of valuable liquids to designated analysis centers on a solid surface as part of 
microscopic chemical factories 0]. Besides the numerous experimental challenges associated 
with these systems there is also the theoretical challenge to understand the corresponding 
highly inhomogeneous fiuid structures (and ultimately the fiow dynamics) as well as to guide 
the design on the basis of this insight [^. 

In view of this goal the study of the fiuid structures in a single wedgelike groove with 
opening angle a serves as a paradigmatic first step. For a = vr the geometry reduces to 
the well studied case of a planar substrate. For decreasing values of a the fiuid is squeezed 
whereas for a > vr a ridgelike solid perturbation is projected into the fiuid. Due to the 
unbound accessible space for all values of a the bulk properties of the fiuid are unchanged 
by the presence of the wedge. (This is an important difference as compared with pores 
in which the confinement alters the bulk properties.) At low temperatures one finds pro- 
nounced packing effects for the fiuid particles near the corner of the wedge which differ 
significantly from those encountered at planar surfaces p|-pn|. Even stronger effects occur 
if the wedge is exposed to a vapor phase. In this case a liquid-like meniscus is formed at 
the bottom of the wedge which undergoes a filling transition at a temperature Tq, below the 
wetting transition temperature Ty/ of the corresponding planar substrate ||T^-[T5|. It has 
been demonstrated that X-ray scattering experiments at grazing incidence are capable of 
resolving such interfacial structures |lT6| , p!7| . 

If the temperature is increased sufficiently above the triple point of the fiuid one en- 
counters its bulk critical point T^. This can be either the liquid-vapor critical point of the 
fiuid or the demixing critical point if the fiuid is a binary liquid mixture. The experience 
with planar surfaces tells |]TB|,|T^ that at bulk criticality the confinement triggers interesting 
new surface critical phenomena. In this case the local order parameter is perturbed near 
the surface within a layer whose thickness is governed by the diverging bulk correlation 
length ^±{t ^ 0) = ^o'l'^l where t = {T — T^/Tc] v denotes the universal bulk critical 
exponent and is the nonuniversal amplitude above (+) and below (— ) T^. Thus at the 
perturbation due to the surface intrudes deeply into the bulk. 

In order to extend our knowledge of the structure of confined fiuids to these elevated 
temperatures we set out to investigate the aforementioned surface critical phenomena in a 
wedge. Based on analytic calculations pOhpTj and computer simulations |p8| , ^ the local 
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critical behavior in a wedge has been analyzed for the case that the corresponding planar 
surface exhibits the so-called ordinary or special surface phase transitions ||l^,l^. In the 
corresponding magnetic language the ordinary phase transition corresponds to the case that 
the couplings between the surface spins remain below the threshold value of the multicritical 
special transition beyond which the surface can support long-ranged order even above Tc 
||l8| , p^ and in which there are no surface fields. The ordinary transition has also been studied 
for different shapes of the confinement such as parabolic ones ppj-p^]. The main concern of 
these studies has been critical edge exponents describing the leading critical behavior of the 
order parameter near the corner where it differs from that at planar surfaces and in the bulk. 
One finds that in general the edge exponents depend on the opening angle a. These studies 
aim at describing magnetic systems in the absence of surface fields as they are experimental 
accessible, e.g., by scanning electron microscopy p3[ . 

However, the aforementioned surface and edge universality classes are not applicable 
in the present context of confined fluids. For example, in a binary liquid mixture one of 
the two species will prefer the confining substrate more than the other which results in an 
effective surface field acting on the order parameter given by the concentration difference 
between the two species. Similar arguments hold for a one-component fluid near liquid- 
vapor coexistence. These surface fields give rise to another surface universality class, the 



so-called normal transition |18,19|. It differs from the ordinary and special transitions 



discussed above in that these symmetry breaking surface fields generate a nontrivial order 
parameter profile even if the bulk is in the disordered phase, i.e., for T > Tc. For a planar 
surface this gives rise to the well studied so-called critical adsorption phenomenon (see, e.g.. 



Refs. [p!8| , p!9| , p4| -p^.) Here we extend the corresponding field-theoretical analysis for the 
normal transition to the wedge geometry under consideration. After presenting our model 
and discussing general scaling properties (Sec. y) we analyze the order parameter profile in 
Sec. m and the structure factor in Sec. By focusing on the critical temperature T = Tc 
we are able to obtain analytical results. They are summarized in Sec.jV]. Appendices ^and 
contain important technical details needed in Sees. |TT| and 

Experimentally, wedges are characterized by a finite depth A. Moreover, in many cases 
they are manufactured as a laterally periodic array of parallel grooves and ridges. The 
structural properties of a fluid exposed to such a periodic array sensitively depends on the 
ratio A/^ of the depth and the correlation length. For A/^ <^ 1 the system will resemble 
a critical fluid exposed to a mildly corrugated substrate which is expected to give rise to 
corrections to the leading critical adsorption behavior at a planar substrate. For A/^ ^ 1 
near the tip of a ridge and near the bottom of a groove the presence of their periodic 
duplicates becomes irrelevant so that in this limit the edge singularity of a single edge with 
a > n and a < tt, respectively, will prevail. Our present study deals with the limit A —>■ oo 
and ^ — >■ cxD taken such that A/^ ^ 1. Certainly, in a later stage it will be rather rewarding 
to study also the crossover regime A/,^ ^ 1 in which the grooves and the ridges as well as 
the way they are joined together play an equally important role. The occurrence of the 
unbending transition in the case of wetting on corrugated substrates indicates that the 
regime A/^ ^ 1 might exhibit rather interesting new phenomena. 
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II. MODEL AND GENERAL RESULTS 



Our investigation is based on the Ginzburg-Landau Hamiltonian Ti, determining the 
statistical weight exp(— 7i) for a scalar order parameter which represents a critical system 
within the Ising universality class. According to the wedge geometry, 7i is given as a sum of 
three contributions, i.e., H = Hb + Hs + 'He, which refer to the bulk (b), the surface (s), and 
the edge contribution (e), respectively. In cylindrical coordinates x = (r, 6',X||) (see Fig.|l]) 
the individual contributions are given by 

nb[<P] = lyej^^drr I d'^-2^|||i(V0)^ + ^0^ + |0^} , (2.1a) 



ns[<P\ = j'^ dr j d'^-2a;||{|[0(r,O,X||)]2-/ii0(r,O,X|,) 

+ ^[0(^">X||)]2 - /i2 0(r,a,X|j)| 



(2.1b) 



j rf"-'xi|||[0(O,O,X||)]2-/i,0(O,O,X||)} , (2.1c) 

where d is the space dimension and xy parameterizes the (d— 2)-dimensional subspace parallel 
to the edge along which the system is translationally invariant. The bulk Hamiltonian H.^ 
given by Eq. ( |2.1a| ) represents the standard Ginzburg-Landau 0^ model in the absence of 
external bulk fields. The bulk parameter r is proportional to the reduced temperature 
t = (T — Tc)/Tc, g is the bulk coupling constant, and a denotes the opening angle. The 
surface Hamiltonian Tis given by Eq. (|2.1bD captures the effect of the two semi-infinite planar 



surfaces forming the geometric boundaries of the wedge located at ^ = and 6 = a, 
respectively. Note that Tib and Tis take their standard fixed-point form with a surface 
enhancement c,- and a surface field hj for each of the surfaces i = 1,2. Cubic surface fields 



3J] will be disregarded here. The surfaces meet at the opening angle a along the edge of 
the wedge which gives rise to the third contribution Tie- The edge Hamiltonian Tig given 
by Eq. ( p.lc|) has the same structure as the surface contribution Hs and is characterized by 
an edge enhancement Ce and an edge field fl^. The total Hamiltonian Ti constitutes a 
renormalizable model (see Subsec. IVA.3 in Ref. [p!9| ). 

In this investigation we are exclusively concerned with the normal transition which is 
characterized by nonzero values of all surface and edge fields hi,h2, and he and arbitrary 
values of the surface enhancements Ci and C2 and of the edge enhancement Cg. For nonzero 
fields hi or /i2 the system is ordered at any finite point even for T > T^. Asymptotically close 
to the critical point the universal properties of the corresponding order parameter profile 
are linked to the critical adsorption fixed-point of the corresponding renormalization group 
description. The ensuing scaling functions refer to the scaling limit r —>■ oo and ^ — > cxo, 
where the ratio r/^ is kept fixed forming a finite scaling variable. At the critical adsorption 
fixed-point the surface fields are infinitely large so that the order parameter profile 

M{r,9,t-a) = (0(r,0,X||)) (2.2) 
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diverges at the surfaces 9 = and 9 = a of the wedge according to a power law. (We recall 
that such divergences refer to the renormalization group fixed-point whereas actually the 
divergence of the order parameter profile is cut off at atomic distances from the surfaces.) 
Consequently, at this so-called normal transition the surfaces can differ at most with respect 
to the sign of the fields. In the following we assume that the surface fields have the same 
sign and we denote this configuration as (+, +). Note that only Hb in conjunction with the 
aforementioned boundary conditions is left to determine the functional form of the order 
parameter profile. 

Close to the critical point the order parameter profile takes the scaling form 

M{r,9,t;a) = a\tfP±{r/^±,9;a) , (2.3) 

where ^± = ^o'l^T'^ with t = {T — Tc)/Tc ^ is the correlation length and -P±(C±5 9] a) are 
the corresponding scaling functions with the scaling variable 

C± = r/U ■ (2.4) 

The amplitude a in Eq. (|2.3| ) is the nonuniversal amplitude of the bulk order parameter 
Mf, = a|t|^, T < Tc, where (3 is the corresponding bulk critical exponent. The scaling 
functions P± are universal and have the limiting behavior P+{C+ oo,9]a) — >■ and 
P-{C- oo,9]a) — > 1, respectively. Note that P+ and P_ are universal but depend on 
the definition of the correlation length because ^± enters the scaling argument (±. In 
the opposite limit C± 0, i.e., T —>■ T^, the scaling functions P± exhibit short- distance 
singularities in form of power laws which reflect the anomalous scaling dimension of the 
order parameter 0: 

P±{C± ^ 0, 9; a) = C±i9; a) Cj'"" . (2.5) 

The ratio /5/z/ of critical exponents has the values 1 in (i = 4, ^ 0.5168 in ci = 3 [4C], and 
1/8 in (i = 2. Equation ( ^.51 ) implies a power law dependence on r for the order parameter 
profile at criticality: 

M(r,e,t = 0;a) = aC^{9;a){r H^Y^I'' . (2.6) 

The scaling functions C±{9;a) appearing in Eqs. ( |2.5| ) and (|2.6| ) are universal but depend 
on the definition of the correlation length because they are derived from P±. (However, the 
product C±{^^)^^'^ is invariant with respect to different choices for the definition of the bulk 
correlation length.) In order to be specific we choose as the definition for ^ the so-called 
true correlation length which governs the decay of the two-point correlation function in the 
bulk system |3^. Equation (p.6|) implies the relation 



C^i9;a)/C^i9;a) = {^^/^^r^/'' (2.7) 

between the scaling functions C± involving the universal ratio Co'/Co'- 

In the limit a = n the wedge geometry coincides with the semi-infinite geometry. In this 
limit the universal scaling functions C± reduce to 

C±{9;a = 7r) =c±{sm9)-^/'' (2.8) 
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where the universal amphtudes c± with c+/c_ = (■Co'/^o ) govern the critical adsorption 
profile at a planar surface |^ so that 

M(r, 0, t = 0; a = tt) = M^i2{z = r sin ^, t = 0) = a c± {z/i^yf^'" (2.9) 

yields the order parameter profile M00/2 in the semi-infinite system at criticality. The profile 
Moo/2 (-^0) t = 0) at a reference distance zq from the surface can be used in order to construct 
the ratio 

M{r,d,t = Q]a) , 

which is independent of the nonuniversal amplitudes a and Here we have introduced the 
length ratio p = v/zq and the scaling function 

C{e;a) = C+{e;a)/c+ = C^{e;a)/c^ . (2.11) 

We note that C is not only universal but also independent of the definition for the correlation 
length because the Ihs of Eq. (p.lO|) and p are expressed in terms of order parameter profiles 
and distances, respectively, without resorting to the notion of the correlation length. 

In the limit that the normal distance z = r sin 6 from the surface of the wedge is much 
smaller than r, i.e., for ^ ^ with r fixed, the order parameter profile reduces to the planar 
semi-infinite behavior 

M{r,9,t = 0;a) = a c± (z/^^)"^/'^ ^ a c± r^/'^(r/eo^)-^/'^ , 9^0 , (2.12) 

so that Eqs. (^) and j ^^AUj ) imply 

C{e ^0;a) ^ 0-^/^ , C{e ^ a;a) ^ {a-e)-^''' , (2.13) 

where we have used the symmetry property 

C{e-a) = C{a-e;a) . (2.14) 



III. ORDER PARAMETER PROFILE AT Z 



c 



At the critical point the radial dependence of the order parameter profile in the wedge is 
given by the power law r~^/'^ [see Eq. ( ^.6| )]. The dependence of the profile on 9 is captured 
by the universal amplitude function C{9;a) [see Eq. (|2.10| )], which in c? = 2 is determined 
by conformal invariance arguments (see Sec. p.lllJ| ). In > 2, however, one has to resort 



to explicit field-theoretical calculations, assuming that the amplitude function C, as other 
universal quantities, is a smooth function of the spatial dimension d. In the following we 
focus on the corresponding mean-field description, i.e., lowest order perturbation theory 
within the field-theoretical approach. The mean-field results become exact at the upper 
critical dimension dyic = 4 of the field theory described by Eq. ( p.l| ), i.e., for d y 4. In order 
to provide an estimate for the order parameter profile in d = 3 we shall use the exact results 
in d = 2 and the mean-field results in c? = 4 for an interpolation scheme between d = 2 and 
d = 4, where the correct scaling arguments and short- distance singularities are implemented. 
For the order parameter profile studied here we shall use the value /3/z/ ~ 0.5168 [^ in = 3 
instead of the mean-field value P/u = 1. Thus only the amplitude function is treated in 
lowest order perturbation theory. 
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A. Mean-field theory 



Within mean-field theory it is convenient to introduce the reduced profile 

m{r,9,T;a) = ^/gju M{r,6,T;a) ; 
m(r, 9, r; a) is determined by the Euler-Lagrange equation 



92 



+ 



Id 1 92 

+ 



m 



rm + 2 m 



(3.1) 



(3.2) 



for r > and < ^ < a. The divergences of the profile m near the surfaces and the edge 
of the wedge account for the effect of the surface and edge parts of the Hamiltonian in the 
(+, +) configuration considered here. For r = Eq. ( |3.2| ) implies for the universal amplitude 
function 



C{9;a) = rm{r,9,T = 0;a) 
[see Eq. ( p.lO|) , here P/u = 1] the equation 



d9^ 



C{9-a) = -C{9;a) + 2[C{9;a)Y 



(3.3) 



(3.4) 



Equation ( p.4|) leads to exact results for C in the limit d y A. Note that Eq. ( p.4|) can 
readily be solved using the order parameter profile obtained for a film geometry with (+, +) 
boundary conditions Specifically, Eq. ( |3.4| ) is equivalent to the Euler-Lagrange equation 
for a film with thickness L as given by Eq. (A3) in Ref. if the rescaled temperature rL^ 
there is replaced by — a^. The expression 9/a for the wedge corresponds to z/L for the film. 
The solution of Eq. (|3.4| ) can be expressed by the Jacobian elliptic functions sn and dn, 



C{9;a) 



1 dn(^/v/l -2F(a);fc(a)) 

v/1 - 2k^{a) sn{9/y^T^^W{a);k{a)) 



< a < TT 



(3.5a) 



in terms of the modulus k of the complete elliptic integral K{k) of the first kind (see 
Appendix ^) where k = k{a) is implicitly given by 



a = 2K{k)Vl - 2A;2 , 1/2 > A;^ > , 
(see Eqs. (A12) and (A13) in Ref. g) and 

C(9]a) = , ^ , ^ , 7i<a<2n 

v/1 + A;2(«) sn{9/^l + k^{ay,k{a)) ~ ~ 

where now k = k{a) is implicitly given by 

a = 2K{k)Vl + k^ , < A; < k^ax ^ 0.90953 , 



(3.5b) 



(3.6a) 



(3.6b) 



(see Eqs. (A14) and (A15) in Ref. |^). In the following we consider only the case < a < vr; 
according to Eq. ( p.6|) the extension to the case vr < a < 27r is straightforward. Note that the 
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parameterizations of the opening angle a as given by Eqs. ( p.5b|) and (|3.6b| ) are monotonous 
functions of the modulus k and therefore the inverse function k{a) is uniquely defined but 
cannot be expressed in closed form in general. The special case a = vr corresponds to 
k{a. = tt) = so that C{Q\ a = vr) = 1/ sin^, which reproduces the mean-field result for the 
order parameter profile in the semi-infinite geometry at the normal transition expressed in 
polar coordinates [see Eqs. ( |2.8|) and (|2.11| ), here (3/v = 1]. The shape of C{9\ a) for d = 2 
[see Eq. ( |3.24| ) in Sec. [lllL)[ | and d = 4 [see Eq. (|3.5|)] for fixed opening angles a = n/2 and 
a = 7r/4 is shown in Fig. |^ where the divergences at 6' = and 6 = a according to Eq. ( p.l3| ) 
have been split off. As an estimate for ci = 3 we show in Fig. ^ the linear interpolation 
between the corresponding curves for d = 2 and d = 4. 



B. Contour lines of the order parameter profile 

One of our motivations for the investigation of criticality in the wedge geometry is the 
frequent use of grooved substrates as templates to impose inhomogeneous structures on 
liquids. In the wedge geometry the density (or concentration) of the adsorbed critical fluid 
is inhomogeneous within planes perpendicular to the edge. The intersection of these planes 
with surfaces of constant order parameter M(r, 6,t = 0;a) = Mq renders the contour lines 
r = r{6] a, Mq) of the order parameter profile, which characterize the inhomogeneous shape 
of the order parameter configuration in the wedge. Far away from the center of the wedge the 
order parameter profile approaches the profile for a planar substrate so that the contour lines 
asymptotically become straight lines parallel to the surfaces forming the wedge (compare 
Fig.|l]). As expected, the deviation of the contour lines from these parallel straight lines is 
largest in the center of the wedge. The normalized order parameter profile defined by the 
Ihs of Eq. (|2.1CI|) yields the contour lines 



ma)r' , (3.7) 



where for a given contour line the reference distance zq is fixed by the asymptotic distance 
of this contour line from the wedge surface, i.e., by Mq = Moo/2{zo,t = 0) (see Fig. p. 
Thus the angular dependence of every contour line exhibits the same universal shape in 
units of the corresponding asymptotic distance Zq. Equations (|2.13|) and (|3.7|) imply that 
r{9 — > 0, a) ~ independent of d and a. 

Within mean-field theory (for which f3/u = 1) we obtain from Eq. ( p.5|) in the center of 
the wedge, i.e., for = a/2: 



r{e = a/2,a) VI - k'ja) 

= C{e = a/2;a) = d = A. (3.8) 

Zq y/l — 2k''{a) 

According to Fig.|^ the distance R of the intersection of the asymptotes from the edge is 
given by R = Zq/ sin(a/2) which implies for the deviation Ar = r — i? of the contour line at 
the center of the wedge 



Ar 



v/1 - k^{a) _ 1 
Vl - 2k^{a) sin(a/2) 



zo , d = 4: . (3.9) 
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In order to discuss Eq. ( p.9|) we first consider a very wide wedge, i.e., the hmit a — > vr for 
which k'^ — > 0. In this limit the deviation Ar decays according to [see Eqs. ( p.5b|) and (A2)] 



In the opposite limit of a very narrow wedge, i.e., for a — > 0, the deviation Ar of the contour 
line in the center diverges according to 

Ar O 

— = - + 0(a) , d = 4 , (3.11) 

zq a 

where = \/2K{l/\/2) — 2 ~ 0.622. This shows that the contour lines are expelled from 
the wedge into the bulk for small opening angles. The overall behavior of Ar for < a < vr 
for d = 2 [see Eq. (|3.26|) in Sec. |11ID| and d = A [see Eq. ( p.9|) ] is shown in Fig.|^. Again, 



the curve for d = ?> represents the estimate obtained by the linear interpolation between the 
curves for d = 2 and d = A (compare Fig.ET). 



C. Distant wall corrections 



Far away from the center of the wedge the contour lines of the order parameter profile 
approach straight lines which are parallel to the wedge surfaces (compare Fig. p. This 
implies that in this limit the deviation of the order parameter profile from the corresponding 
profile in a semi-infinite geometry becomes small. This small deviation can be interpreted 
in terms of the so-called distant-wall correction which is obtained from the short- distance 
expansion [P^ -|46[| of the order parameter near, say, the surface ^ = of the wedge. For any 
scaling operator \E' the short- distance expansion at the normal transition can be written as 



^^(r,^,X|| 



(^(r,^^,X||))„o/2 [l + 6*ree(r,0,X||)(r^)^ + 0(6^+^)] 



(3.12) 



where Tgg = T±± denotes the stress tensor component perpendicular to the wall at ^ = 
and r6 <^ ra is the small distance from the wall. For the order parameter profile M = {(p) 
we infer from Eq. ( p. 121 ) 



M{r,e;a) = M{r,e;7r) [1 + 6'' + ...' 



(3.13) 



in analogy to the short-distance expansion of the order parameter profile in the film geometry 
p7| , where = b^{T0g)r'^. From the general expression of the stress tensor at criticality 
(see Eq. (3.1) in Ref. |HB|) we obtain in polar coordinates 



To, 



2rAde 2\dr 2 ^ " 



2 g^ d-2. 



(3.14) 



up to contributions of the order g^, where Vy and Ay denote the components of the gradient 
and the Laplacian along the edge, respectively. 

Within mean-field theory the thermal average of Eq. (|3.14|) is obtained by setting d = A 
and replacing by {(f)) = ^JvIJg m [see Eqs. (|3.1| ) and (p. 31)]. Since the result is independent 
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of 9 we can use the fact that the order parameter profile is symmetric around the midplane 
6' = a/2 in order to obtain 



(T,,) = ^ [C{a/2;a)r-[C{a/2-a)Y 



d 



From Eq. p.5| ) we thus find 



{Tee) = - 



6 e{a)[l-k\a)] 
^ [l-2P(a)]2 



d 



(3.15) 



(3.16) 



which exactly corresponds to (Tzz) in the film geometry with (+, +) boundary conditions 
(see the first line of Eq. (3.4) in Ref. g for y = -I). If Eq. j ^J^ is inserted into Eq. ( p32D 
for = (f) the structure of the distant-wall correction as given by Eq. (|3.13| ) is confirmed. 
However, for a complete identification the amplitude b^p must be determined. Within the 
present mean-field approach Eq. ( |3.13|) can directly be verified by expanding Eq. ( |3.5| ) for 
— > with a fixed. From the expansions of the Jacobian elliptic functions [see Eq. ([^)] we 
obtain 



m(r, 9] a) = m{r, 9; vr) 



1 + 



[1 -2fc2(a)]2 



9^ 
10 



105 



+ 0{9^) 



, d = 4, (3.17) 



from which one can read off = —g/60 within mean-field theory. The next-to- leading 
correction 0{9^) in Eq. ( p.l7| ) is generated by the operator ArTee where denotes the 
radial part of the Laplacian in polar coordinates (r, ^,X||). In general the corresponding 
correction is 0{9'^^'^) as indicated in Eq. (|3.12|) . In principle AyTee also yields a contribution 
to the short- distance expansion, but one has Ay (Tee) = due to the translational invariance 
along the {d — 2)-dimensional edge. We note that in the film geometry the counterpart 
of the operator ArTee in the wedge does not contribute to the next-to-leading distant-wall 
correction of the order parameter profile due to the translational invariance in the film. 
Therefore in the film the next-to- leading correction is given by the operator {T±±y, which 
gives rise to 0{9'^'^) corrections only. The latter operator appears also in the wedge geometry 
and gives rise to the 0{9^) corrections indicated in Eq. (|3.17|) . 

The distant- wall correction has an interesting geometric interpretation. The contour lines 
of the order parameter profile in the wedge and in the semi-infinite geometry, respectively, 
belonging to the same value Mq of the order parameter are displaced by an amount S{y; a) 
where S + zq = z = r sin 9 and the distance y is related to 9 by (compare Fig. |I|) 



y 



y{9,a) = r{9, a) cos 9 - [R{a) + Ar{a)]cos{a/2) . 



Within mean-field theory one obtains for the above choice of Mq [see Eq. 
and using t/zq = C{9; a) 



5{y{9,a);a) 



C{9] a) sin(0) 



m(r, 9] a) 
m{r, 9; vr) 



- 1 



d = A 



(3.18) 
here l3/u = 1] 



(3.19) 



According to Eqs. ( |3.13D and ( |3.19| ) the ratio 6{y{9, a); a)/zo as function of 9 represents the 
distant-wall correction of the order parameter profile. Its leading terms can be read off from 
Eq. ( |3.17|) . For general d one has [see Eq. (|3.7| )] 
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[C{e,a)Y/^ sm{e) - 1 



(3.20) 



so that 6/zo (l/6)[(7r/a)2 - 1] 9^ for = 2 and ^ ^ [see Eq. d^M) in Sec. pTD|] . The 
behavior of 6{y, a) as function of y for fixed opening angles a = 7r/2 and a = n/A is shown 
in Fig. H for d = 2 [see Eq. (|3.24 )] and d = 4 [see Eq. ( |3.19 )]. In order to obtain an estimate 
for (i = 3 we introduce 



zo \a J 



(3.21) 



e 



-1 



One has 6/zq ~ 6'^ for 6 —>■ or ^ y^'^ for y oo due to y{6 ^ 0, a) = r{6 —>■ 0,a 
[see Eq. ( p.l8| )]. Therefore f{9;a) tends for ^ to a constant value (depending on d). 
For y = one has 29/a = 1 and / = Ar sin(a/2)/2;o (compare Figs.|T] and |^). A reasonable 
estimate for S{y; a) in c? = 3 can be obtained by interpolating f{6; a) linearly between d = 2 
and d = 4 and using this approximation for / in 



S{y{e,a);a) 



'20 

a 



d = 3 . 



(3.22) 



The curves for c? = 3 in Fig.^ correspond to the rhs of Eq. ( p.22|) for fixed values a = 7r/2 
and a = n/A. For y = one has 26 /a = 1 so that for y = the linear interpolation of / 
implies a linear interpolation of 6 [see Eq. ( |3.21| )]. 



D. Exact results in d 



At criticality systems exhibit not only scale invariance but, more generally, conformal 
invariance This property is particularly useful in d = 2, where the large conformal group 
provides mappings between many different geometries MM . In higher spatial dimensions only 
Mobius transformations are available as conformal mappings (small conformal group) which 
map geometries bounded by planes and spheres onto other geometries bounded by planes 
and spheres. In c? = 2 a wedge with opening angle a can be conformally mapped onto a 



half-plane, for which the profile is known. This leads to [^ 



M{r,e,t = 0;«) = {(t))^edge = A (Tc/af^" [rsmiOn/a)]-"/" 



(3.23) 



in the wedge, where P/u = 1/8 within the Ising universality class considered here. Note 
that conformal invariance completely determines the universal amplitude function C{6; a) 
[see Eq. ( CT) ] as 



C{9; a) = (n/af/" [sin(^7r/a)] 



d 



(3.24) 



which would remain unspecified by scale invariance considerations alone. 
Using Eqs. ( p.24|) and ( p.7|) we obtain the relation 



Zo 



Ti/a 
sin(6'7r/a) 



d 



(3.25) 
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for the contour lines of the order parameter profile. The deviation Ar from the asymptotic 
contour lines at the center of the wedge, i.e., 9 = a/2, is defined as shown in Fig.|l|. We 
obtain 

Ar 7r 1 



7r 
a 



zq a sin(a/2) ' 
Eq. {^^) yields the expansion 

-2. . 2 



d = 2 . 



1 



TT 



1 - 



a 



TT 



1 



d 



In the limit a 
Ar 

Zq 

which is the analogue of Eq. ( |3.1CI|) for c/ = 2. For a we obtain from Eq. ( p.26| ) 

Ar 7r-2 , 
— = + 0{a) , d = 2 , 



TT 



(3.26) 



(3.27) 



(3.28) 



zq a 

which is the analogue of Eq. ( PH) . Note that Eqs. (|3l0|) and ( pTj) and Eqs. (PD and 
), respectively, are very similar even in a quantitative sense, i.e., the dependence on the 



spatial dimension d is weak. 

We close this subsection with a discussion of the short- distance expansion of the order 
parameter according to Eq. (|3.12| ) for the case d = 2. The thermal average of the stress tensor 



(T) in the wedge can be obtained by means of a conformal mapping from the half-plane |25 
and reads in polar coordinates 

1 



{Trr) — ~{Tee) 



247T 



/TT\ 2 

Vo;/ 



1 



(3.29a) 



{Tre) = {Tor) = , (3.29b) 

where c is the conformal anomaly number. The distant-wall correction for the order parame- 
ter profile according to Eq. ( |3.13| ) can be directly calculated from the expansion of Eq. ( p. 23 ) 
for small 6 with the result 

M{r,e,t = 0;a) = M{r,e,t 



0:tt) 1 + 



- 1 



(3.30) 



From Eqs. ( pJ2D , (1091 ), and ( p30|) we obtain = -{Att/c) f3/u in d = 2. For the Ising 
universality class considered here (c = 1/2 and /5/z/ = 1/8) we have = —tt. Finally, 
we note that the next-to-leading distant-wall correction 0{6'^) indicated in Eq. ( |3.3U| ) stems 
from ArTge and {TegY because both operators have the same scaling dimension four. 



IV. STATIC STRUCTURE FACTOR 

The static structure factor 5'(|q|) in a bulk system is given by the Fourier transform of 
the two-point correlation function Gdr — r'|), where q denotes the momentum transfer and 
r and r' are positions in space. Due to translational invariance and isotropy of the bulk the 
moduli of the momentum transfer and the distance in space remain as the only arguments 
in Fourier space and real space, respectively. In the wedge geometry translational invariance 
and isotropy only hold in the (d— 2)-dimensional subspace along the edge so that the position 
dependence of G is more complicated. Apart from modifications of the eigenmode spectrum 
caused by different boundary conditions we follow the derivation of G for a wedge as given 
inRef. 
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A. Mean-field theory 



Starting from the formal definition ^(r; r') = [6m{r) / 6h{r')]h=Q , wliere m(r) denotes the 
order parameter profile and h{r) is an external spatially varying field, we obtain 



- AG(r; r') + Qm\r)G{r; r') = 6{r - r') 



(4.1) 



at the critical point and within mean- field theory. For r = (r, ^^,X||) we apply a Fourier 
transform with respect to xy and define ^(p; r, 6; r', 9') by 



G(r,^,X||;r',^',x;i 



(27r)'^-2 



5(p;r,^;r',^')exp[ip(x|| -xj,)] . 



(4.2) 



From Eqs. ( ^.1| ) and ( ^.21 ) and by using m(r) = C{9\ a)r ^ [see Eq. ( ^.3| )] one obtains 



Id I 

+ -7^ + 



^ ^2 QQ2 



V 



(4.3) 



We solve Eq. (|4.3| ) in terms of eigenf unctions to a spectrum E of eigenvalues for the operator 
on the left hand side of Eq. ( [4.3| ). Each eigenfunction can be written as a product of a radial 
part R{r) and an angular part ip{0). For the shifted spectrum = E — p"^ the corresponding 
eigenvalue equations are given by 



and 



<'(r) + r-i<(r) + {k^ - \nr~^)Rn{r) = 



(4.4) 



(4.5) 



where Eq. j WA ) is a Bessel equation with parameter k and Eq. (|4.5|) is a Lame equation. 
In general the discrete eigenvalues cannot be given in closed form (see Appendix ^ for 
details). Apart from the specific spectrum {A„} and the corresponding eigenfunctions {V'n} 
one can represent the solution of Eq. ( [4.3| ) in the same form as in Ref. pO 



I.e., 



2 °° r-oo 

S'(p; r, 9; r', 9') = — / dn k 



J.n-(nr)J, 



Kr 



n=3 

oo 



+ p^ 



-Y.^V>rSP^<)Kv>rSP^>) ^n(^)^n(^^') 



(4.6) 



?i=3 



where p = |p|, r< = min(r, r'), and r> = max(r, r'). In order to determine the decay of the 
two-point correlation function in real space we insert Eq. ( [4.6| ) into Eq. ( [4.2| ) and proceed 
along the lines described in Ref. |2^. The decay of G away from the edge is governed by 
the critical exponent rig±{a) according to 



G{r, 9, X|| ; r' ^ oo, 9' , Xn 



(^^/^d-2+r;ex(a) 



(4.7) 
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and as in Ref. [^] we find from Eq. ( [4.2| ) tliat in the limit r'/r oo the smallest eigenvalue 
in the spectrum {A„} determines the exponent rie±{a). With the notation used in Appendix 
B the smallest eigenvalue is A3 and in accordance with Ref. we find r]e±{a) = y/Xs- 
Along the edge, i.e., in the limit |x|| — xjj| — > 00 for r, r' > the decay of G is governed by 
the critical exponent ?7e||(«) according to 

G(r,^,X||;r',0',xj, ^00) ~ ^^^^ _ ^,^,-2^,,,^.^ ■ (4-8) 

In agreement with the general scaling relation 

^ex(«) = ^^^f^ , 2<d<A, (4.9) 

we find rie\\{a) = 2rie±{a) within mean- field theory (for which r] = 0). From the approximate 
evaluation of the spectrum {A„} [see Eq. (p6| ) in Appendix B] we obtain 




Ve±{a) = ^ -^/l-^C(«/2) + ^(7r2 + ^a2_2e(a/2)) (4.10) 

where ({u) is the Weierstrass ^-function [see Eq. (|A9| )] and 

, , Ak^(a)\l- k'^(a)] 4 

= |l-2P(a)p + 3 ("^' 

[see Eq. (PT|)]. In the limit a n corresponding to a planar surface Eq. ( ^4.10| ) becomes 
exact and one obtains 7]j_ = i]f.±{TT) = 3 in accordance with the mean- field result for the 
semi-infinite geometry at the normal transition |TH]. The dependence of rif.± on a is displayed 



in Fig.^ for d = 2 [see Eq. (|4.16| ) in Sec. [IVB[| and d = 4 [see Eq. ( [4.10|) j, where the overall 



1/a dependence has been split off. The estimate for (i = 3 is obtained by interpolating 
linearly between the curves for d = 2 and d = 4 (compare Figs.^ and |^). The exponent 
rie±{a) increases monotonically as a decreases. This implies that at the normal transition 
the correlation function decays more rapidly away from an edge {a < vr) than away from a 
planar surface (a = vr). 

Finally, we note that from the general scaling relation for the edge exponent 

^ , , uld — 2 + riew (a)] 

/5e(a) = ^ 2<rf<4, (4.12) 

which describes the singular temperature dependence of the order parameter at the edge, 
one has within mean- field theory [see Eq. ([4.9|) , here 77 = and u = 1/2] 

A(„) = i±MM, , = 4. (4.13) 

In the limit a — > tt corresponding to a planar surface Eq. ( [4.13| ) yields /5e(7r) = jSi = 2, in 
accordance with the general scaling relation jSi = dv for 2 < d < 4 ai the normal transition 
for the semi-infinite geometry p6|. 
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B. Exact results in d = 2 



In d = 2 rigorous results for the edge exponent i]e±{a) can be obtained starting from the 



correlation function in the half-plane at the normal transition p9 

' {x - x'f + y'^ + y''^ 



G{x,y]x',y') = {yy') 



yy 



(4.14) 



where x, x' are the coordinates of the two points along the surface and y, y' are their distances 
from the surface |^9[. Applying the conformal mapping and carrying out the limit r'/r oo 
at fixed r > leads to 



G(r,^;r',e') 



7r\^ 
a 



TT 



a 



sin f —9 ) sin f —9' 



TT 



a 



(^ir'?)/2 



rr 



V2 (L 



-7r7y||/(2a) 



(4.15) 



for the correlation function in the wedge. From Eq. ( [4.15| ) one can read off the scaling 



relation 49 



'7e±(a) 



2 a 2 



(4.16) 



In the limit a — >■ vr corresponding to a planar surface one has ?7e±(7r) = [rj + r]\\)/2 = 7]± 
as expected [see Eq. (^.9|)]. In ci = 2 the edge consists of a single point so that one cannot 
define correlations along the edge. However, if we formally define ?7e|| (o^) by Eq. (^4.9|) we find 
rie\\{<y) = (7r/Q;)?7||. From f3i = 2 we obtain r]\\ = 4 at the normal transition within the 
Ising universality class, which implies [see Eq. (|4.12|) , here z/ = 1] 



/3e(«) 



271 

a 



d 



(4.17) 



for the edge exponent of the magnetization. We note that (3e given by Eq. ( |4.171 ) is four 
times larger than its corresponding value at the ordinary transition [pO| -|2^. 



V. SUMMARY 

We have investigated the universal local properties of the order parameter profile in a 
wedge with opening angle a (see Fig.|I|) for the normal transition. We have obtained the 
following main results: 

1. Near the order parameter is determined by universal scaling functions and the two 
nonuniversal bulk amplitudes a and [see Eq. (^.31)] . At Tc the order parameter profile 
reduces to a power law r~^/'^ in radial direction multiplied by a universal amplitude 
function depending on the polar angle 9 and the opening angle a [see Eqs. ( |2.ti| ) and 
( |2.10| )]. The amplitude function is symmetric around the midplane and diverges as 



9 ^^'^ upon approaching the surfaces forming the wedge [see Eq. (|2.13|) ]. 

We have determined the universal amplitude function C{9, a) within mean-field theory, 
i.e., for space dimension (1 = 4 [see Eqs. ( p.5|) and (|3.6|) ], where the order parameter 
profile in the wedge and for T = Tc can be obtained from the order parameter profile 
in the film geometry for T < [see Eq. (|3.4|) ]. In conjunction with exact results in 
d = 2 [see Eq. ( |3.24|) ] we construct an estimate for C{9, a) for d = 3 (see Fig. 



15 



The amplitude function determines the meniscus-hke contour hues of a constant value 
of the critical order parameter profile [see Eq. (|3.7|) and Fig.|l|. The deviation Ar of 
the contour line relative to its asymptotes (compare Fig.|l]) from the corner of the 
wedge vanishes linearly in the planar limit a — vr [see Eqs. (|3.1CI|) and ( p. 27] )] and 



diverges ~ a ^ for a — > [see Eqs. ( p.llD and (|3.28| )]. Figure ^ presents an estimate 
of the function Ar{a) for d = 3. 

4. The contour lines approach their asymptotes as y~'^ for increasing lateral distance y 
[see Fig.|l| and Eq. ( p.20| )]. This follows from an analysis of distant-wall corrections in 
terms of the stress tensor (see Sec. [Ill C|) . The explicit results for ci = 2 and d = 4 



allow one to construct an estimate for the corresponding behavior m d = 3 (see Fig. 

The decay of the two-point correlation function at Tc away from the edge and parallel 
to the wedge is governed by the critical edge exponents rie±{a) and rj^w (a), respectively 
[see Eqs. ( [4.7|) and ( ^l8|) ]. They fulfill the scaling relation rie±{a) = [r] + rie\\{a)]/2 with 
the bulk exponent rj. Based on the quite accurate estimate in d = 4 [see Eq. ( [4.10[ ) 



and Appendix B] and the exact result in d = 2 [see Eq. ( [4.16| )] Fig. ^ presents an 
estimate for ?7e±(a) in d = 3. Equations ( |4.13| ) and ( [4. IT] ) provide the critical edge 



exponent /?e(«) of the order parameter for the normal transition. Equation ( |4.6| ) gives 
the full structure factor for the wedge geometry and for T = within mean-field 
approximation. 
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APPENDIX A: ELLIPTIC FUNCTIONS 



Here we summarize a few properties of elliptic functions as far as they are needed for 
the derivation of the results obtained in this paper. For further information we refer to the 
literature (see, e.g., Refs. [|5D|-p^). The properties of the Jacobian elliptic functions can 



be derived starting from the Jacobi amplitude am(M; k) which is implicitly defined by the 
incomplete elliptic integral of the first kind: 



u 



am(M; k) 



dip 



k"^ sin Lf 



(Al) 



The complete elliptic integral of the first kind K = K{k) is defined by am(i^'; k) = n/2. For 
the derivation of Eq. (|3.1CI| ) we quote the expansion of K{k) in powers of the modulus k: 



K{k) = I 



(A2) 
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From the first derivative of Eq. ( |Al| ) with respect to u one obtains for k < 1 the relation 

d 



dn(u; k) = — amf-u; k) = \/l — Psn^fw; k) (A3) 
ou 

for the delta amplitude dn(u; k) using the standard notations 

sn('u; k) = sin[am(M; k)] , cn{u; k) = cos[am(M; k)] . (A4) 
Due to am(0; /c) = and am(i^; k) = tx/2 one has 

sn(0; fc) = , cn(0; k) = I , dn(0; k) = I , (A5a) 



sn{K; k) = l , cn{K; k) = , dn{K; k) = Vl - k^ . (A5b) 

The derivation of Eq. ( |3.17| ) is based on the Taylor expansions 

1 + e ^ 1 + Uk^ + k^ , 1 + 135P + 135fc4 + A;6 9 
sn(M; A;) = M — H ?r u' + 0{u^) , (A6a) 

dn(«; A;) = 1 - — + '-u^ ^ '-u^ + 0{u^) . (A6b 

2 4! 6! 

Finally, we quote the relation 

pu ^am(M; k) , 

E{am{u; k), k) = / dn'^{x; k)dx= / -J 1 — k'^ siii^ (p d(p (A7) 

JO io 

between the incomplete elliptic integral of the second kind E{x, k) and the delta amplitude 
dn('u; k). E{k) = E{n/2, k) is the complete elliptic integral of the second kind. 

As demonstrated in Appendix A of Ref. |^T|] the mean-field order parameter profile can 
be obtained from the observation that [C{6] a)]^ is a Weierstrass p-function up to an additive 
constant. The Weierstrass p-function is an elliptic function which is related to squares of 
certain Jacobian elliptic functions |]50|,pl|. It solves the differential equation 



[p'iu)]^ = Ap%u)-g2piu)-gs (A8) 

where g2 and gs are the invariants of p. Note that no term quadratic in p appears on the 
rhs of Eq. (|A^). This condition determines the additive constant in the relation between 
p and which appears as the "potential" in the eigenvalue problem in Eq. ( |4.5| ). For the 
derivation of the spectrum {A„} (see Appendix B) the Weierstrass ^-function is also needed. 
It is the negative integral of p{u) and can be written as 



1 r 
= - - / 

u Jo 



PU 




Jo 





dz . (A9) 

Note that C{u) is not an elliptic function. For the exphcit calculation of the spectrum {A„} 
we finally quote the Laurent series of p{u) and ({u) around u = 0: 



p(u) = 4 + — + — + 0(u^) . (AlOb) 
^ ^ ^ u2 20 28 ^ ' ^ ^ 
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APPENDIX B: EIGENMODE SPECTRUM 



The spectrum of the eigenvalue problem defined by Eq. ([4.5|) can be determined along 
similar lines as in Appendix B of Ref. The Weierstrass function associated with C{6] a) 
[see Eq. (|3.5| )] can be written in the form p{6) = €"^[6] a) — a, where a is a constant to be 



determined. Based on Eq. ( |3.4| ) we aim at obtaining a differential equation for p which is of 
the form given by Eq. (|A8|) . This requirement is fulfilled if a = 1/3 which results in 



2Kdn{2Ke/a;k) 
a sn{2K9/a;k) _ 



1 

3 ' 



(Bla) 



92 



2k 



2\2 



(Bib) 



_ 4fc2(l - fc2) _8_ 
~ 3(1-2P)2 + 27 ■ 

The eigenvalue problem in Eq. ([4.5|) can now be cast into a Lame equation ^2 

- f;{e) + 6 p{e)M0) = ^n(^) 



(Blc) 



(B2) 



where e„ = A„ — 2. As discussed in Appendix B in Ref. the eigenvalue spectrum is given 
by the solution of the two equations 



2a„C(a/2) — a 



p (a^ 



2p(a„) - e„/3 



nin 



(B3a) 



4 - 27^3 
27(72 - 9e; 



(B3b) 



The mode numbers n are integers (see below) and a„ are auxiliary parameters with the 
property a„ ^ for n ^ oo. For large n Eq. ( pHf ) can be solved asymptotically by using 
the expansions quoted in Eq. ( |A10| ). From the expansion for p(a„) we obtain for large n 
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1 + l§«n + 0{al' 
540 " ^ 



which implies the expansion 

2 



C(«/2)a„ -f-a^ + 0{al) 

a a„ loU 



nm 



a 



(B4) 



(B5) 



Equation ( P5|) can be solved with the ansatz a„ = 3a2/(n7r)[l + An ^ + Bn + 0{n ^)]. 
Inserting the solution into Eq. (pi) leads to 
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mvr 



a C(a/2) 



36 



mvr 



[a C{a/2) f - 



+ 0{n- 



(B6) 



An = Cn + 2 yields the desired spectrum. 

The allowed mode numbers n can be obtained by considering the special case a = tt 
corresponding to a planar surface, for which C{9\ vr) = 1/ sin^ and Eq. ( |4.5| ) can be solved 
in closed form. For the present problem the Weierstrass ^-function can directly be derived 
from Eqs. and so that 



a 



dn{2Ke/a; k) 
sn{2Ke/a;k) 



^2 _ 2 2K 

cn{2Ke/a; k) + E{am{2K9/a; k), k) + — d 

3 a 



(B7) 



with K = K{k) as defined in Appendix A. Equation ( [B7D explicitly demonstrates that C(^) 
is not an elliptic function. At the midplane 6 = a/2 Eq. ( p7|) reduces to [see Eq. ( |A5|) ] 



C(«/2) 



2K{k) 



a 



E{k) + 



e-2 



K{k) 



According to Eq. (p.5b| ) the special case a = tx corresponds to = for which ({tt/2) = vr/G 
from Eq. ( P8|) and g2 = 4/3 from Eq. ( pi|) . From Eq. ( p6|) we infer = — 2, i.e., = n^, 
which is indeed the correct eigenvalue spectrum for C{9;tt) = 1/ sin 6*. The eigenf unctions 
are normalizable for n > 3. Therefore A3 is the lowest eigenvalue for this problem. 

For a numerical solution of Eq. ( p3|) using, e.g., the Newton method the asymptotic 
spectrum given by Eq. (P6|) for n > 3 provides excellent initial values for the iteration. In 
fact, these initial values already are within 0.1% of the exact spectrum even for the ground 
state ?T, = 3if0.1<Q;/7r< 1. This implies that the mean-field expression for rie±{a) given 
by Eq. ( [4.10|) is quite accurate if a is not too small. 
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FIGURES 



FIG. 1. Cross section of a wedge perpendicular to its edge with opening angle a. The system 
is translationally invariant in the {d — 2)-dimensional subspace parallel to the edge. The curve 
labeled "M = const" represents a contour line of the order parameter profile (see Sec. IIIB| ). 



r 



and < 9 < a are cylindrical coordinates, z = rsinO is the normal distance from the surface 
of the wedge. The dashed lines are the asymptotes of the contour line M = const which are a 
distance zq apart from the surface. R = zq/ sin(a/2) is the distance between the intersection of 
the asymptotes and the edge of the wedge. Ar + R is the distance between the contour line in the 
center and the edge. The contour line approaches the asymptotes such that 5{y — > oo) — > and 
6{y = 0) = Arsin(a/2) where y measures the distance along the wall. 



FIG. 2. Universal scaling function C{9;a) [see Eq. (|2.10)] as function of the polar angle 9 for 



fixed opening angles a = 7r/2 and a = it/A (compare Fig.[l|). The divergences at = and = a 
are split off [compare Eq. ( |2.13| )]. The curves for d = 3 show the linear interpolation between the 
corresponding exactly known curves for d = 2 [see Eq. (1^] and d = 4 [see Eq. (jsj)]. Note that 



the shape of the curve for d = 2 is independent of a for < a < vr. 

FIG. 3. Universal dependence of the deviation Ar = r — i? in units of the reference distance zq 
of the contour line at the center of the wedge (compare Fig. IT]) on the opening angle a. The curve 



for d = 3 shows the linear interpolation between the exactly known curves for d = 2 [see Eq. (3.26)] 
and d = 4 [see Eq. (p.9|)]. Ar is multiplied by a so that the product aAr attains a finite value for 
a — > [see Eq. ( p. 28]) for d = 2 and Eq. ( |3.11 ) for d = 4; we assume in addition that the power law 



Ar(a ^ 0) ~ a remains valid for 2 < d < 4]. 

FIG. 4. Universal dependence of the displacement 5{y; a) of the contour line from its asymptote 
on the distance y [compare Fig. |^ and see Eqs. (3.20) and ( p.l8| )] for fixed opening angles a = tx 12 



and a = it/ 4. The curves for d = 3 correspond to a suitable interpolation [see Eq. ( 3.22 )] between 
the exactly known curves for d = 2 and d = 4 so that for y/zQ oo they exhibit the correct 
asymptotic decay ^ {y/zo)""^ in accordance with the distant-wall correction. 

FIG. 5. Universal dependence of the edge exponent r?e±(o) [see Eq. (^)] on the opening angle 
a, with the overall 1/a dependence split off. The curve for d = 3 shows the linear interpolation 
between the exactly known curves for d = 2 [see Eq. ( [4.16 )] and d = 4 [see Eq. ( [4.10| )]. The dot 



for a/vr = 1 shows the known value rjx = [t] + ?7||)/2 ~ 2.518 corresponding to a planar surface in 
d = 3. The small difference between the dot and the value of the interpolated curve for d = 3 is a 
measure of the uncertainty associated with the aforementioned linear interpolation scheme. 
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